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ABSTRACT 

The paper is devoted to the methods of determination of the cosmological 
parameters and distortions of ionization history from recent CMB observations. 
We show that the more complex models of kinetics of recombination with a 
few "missing" parameters describing the recombination process provide better 
agreement between the measured and the expected characteristics of the CMB 
anisotropy. In particular, we consider the external sources of the resonance and 
ionizing radiation and the model with the strong clustering of the baryonic com- 
ponent. These factors can constrain the estimates of the cosmological parameters 
usually discussed. We demonstrate also that the measurements of the CMB po- 
larization can improve these estimates and, for the precision expected from the 
PLANCK mission, allow to discriminate a wide class of models. 

Subject headings: theory: - cosmic microwave background - ionization 
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1. Introduction 

Recent observations of the CMB anisotropy such as the BOOMERANG (de Bernardis 
et al. 2000), the MAXIMA-1 (Hanany et al. 2000), the DASI (Halverson et al. 2001), 
the VSA (Watson et al. 2002), the Archeops (Benoit et al. 2002) and especially the 
new CBI data (Mason et al. 2002) and the DASI polarization measurements (Kovac 
et al. 2002) provide the good base for the progress in understanding the most general 
properties of the early Universe. In these experiments, major attention is on the statistical 
properties of the signal and the determination of the power spectrum of the CMB anisotropy 
and polarization. The obtained cosmological parameters for adiabatic cold dark matter 
model (CDM) are well consistent with the Big Bang nucleosynthesis predictions, SN la 
observational data and results obtained from investigations of galaxy surveys. Further 
progress will be achieved with more sensitive experiments such as the launched MAP and 
the upcoming PLANCK missions. They will be able to measure the CMB anisotropy and 
polarization power spectra up to the multipole range I ~ 3 — 4 • 10 3 where the lensing effect, 
Sunyaev-Zel'dovich effect in clusters, secondary anisotropies and some peculiarities of the 
ionization history of the Universe dominate. 

The high precision of the future CMB measurements allows to reveal some distortions 
with respect to the standard model of the hydrogen recombination process occurred at 
redshifts z ~ 10 3 . Here we show that polarization power spectrum contains important 
information about the kinetics of hydrogen recombination and allows to trace the ionization 
history of the cosmic plasma in wide range of redshifts 10 3 > z > 6 — 10. 

For the baryon-dominated Universe the classical theory of the hydrogen recombination 
was developed in Peebles (1968) and Zel'dovich, Kurt and Sunyaev (1968). For the dark 
matter dominated Universe, it was generalized in Zabotin and Naselsky (1982), Jones 
and Wyse (1985), Seager, Sasselov and Scott (1999) (see detailed review in White, Scott 
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and Silk, 1994). Many distortions of the standard model of recombination have also been 
investigated. The delay of recombination due to evaporation of primordial black holes has 
been discussed in Naselsky (1978) and Naselsky and Polnarev (1987). Avelino et al. (2000), 
Battye et al. (2001), Landau et al. (2001) have noted that possible time variations of the 
fundamental constants could also influence the ionization history of the cosmic plasma. 
Sarkar and Cooper (1983), Scott et al. (1991), Ellis et al. (1992), Adams et al. (1998), 
Peebles, Seager and Hu (2000), Doroshkevich and Naselsky (2002) discussed distortions of 
the kinetics of hydrogen recombination caused by decays of hypothetical particles. It is 
worth noting that for such models the energy injection delays the recombination at z ~ 10 3 
and distorts the ionization history of the Universe down to the period of galaxies formation 
at z ~ 5 — 10. 

Recently Naselsky and Novikov (2002) have proposed a model of accelerated 
recombination with non-uniform spatial distribution of baryonic matter. In this model the 
recombination proceeds faster within the clumps and slowly in the intercloud medium. In 
such models at redshifts z ~ 10 3 the mean ionization fraction decreases faster than that in 
the standard model with the same cosmological parameters which mimics the acceleration 
of recombination. Here we show that, for suitable parameters of the clouds, such models 
can improve the best-fit of the observed power spectrum of the CMB anisotropy. 

The non-standard models of hydrogen recombination are characterized by a few 
additional parameters which must be included in the list of the best-fit cosmological 
parameter determination using the CMB anisotropy and polarization data. We call these 
parameters as the "missing parameters" of the theory and show how important these 
parameters may be for the CMB physics. 

In this paper we compare the observed CMB anisotropy power spectrum for the 
reference model of the standard recombination, with models of delayed and accelerated 
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recombination. For the reference model we accept the cosmological parameters from the 
best-fit of the CBI observational data (Mason et al 2002). We show that the models 
considered in this paper can provide better fits for the already observed and the future 
MAP and PLANCK power spectra and can change usual estimates of the cosmological 
parameters. Moreover, the measurements of the polarization allow us to discriminate 
between some of such complex models of the Universe. These results can be important for 
the interpretation of the MAP and upcoming PLANCK data. 

The paper is organized as follows. In section 2 we discuss the generic models of 
distorted recombination. In section 3 we demonstrate distortions of the CMB anisotropy 
power spectrum due to delay or acceleration of recombination. In section 4 we introduce 
the models of accelerated and delayed recombination and discuss the redshift variations of 
the hydrogen ionization fraction in these models. In section 5 we compare the observed 
CMB anisotropics and expected polarization power spectra with ones for the models under 
discussion. Section 6 is devoted to discussion of late reionization and of the polarization 
power spectrum for the delayed recombination models. In Conclusion we make some 
predictions for the MAP and the upcoming PLANCK experiments. 



2. Deviations from the standard model of hydrogen recombination 

The deviations from the standard recombination process caused by the injection of 
additional Ly-a and Ly-c photons at the recombination epoch or by strong small scale 
clustering of the baryonic component can be suitably described in terms of an additional 
source of Ly-a and Ly-c photons which was proposed in Peebles, Seager and Hu (2000). In 
the generalized recombination model the rates of production of resonance, n a , and ionized, 
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rii, photons are described by functions e a (z) and Si(z), defined as follows 

dn a 



= e a (z)H(z)(n b (z)) , (1) 
= e i {z)H(z){n b (z)), (2) 



dt 

dn, 
~di 

where H(z) and (n b (z)) are the Hubble parameter and the mean baryonic density. 



2.1. Delay of recombination 

Peebles, Seager and Hu (2000) considered the models with e a (z) = const, and 
£i(z) = const. However, in general case both parameters e a (z) and Si{z) could depend on 
redshift (see, e.g., Doroshkevich & Naselsky 2002; Naselsky & Novikov 2002). As was shown 
in these papers, for the generation of ionized photons from the decay of super heavy dark 
matter (SHDM) particles we get 

e a (z) ~ £i(z) oc (1 + z)' 1 . (3) 

For models of particle or strings decays the particle number density, decreasing 
function of time defined by the life time of the particles, T x (n x ,t), as follows: 

^ + 3tf(tK = -^— (4) 

dt T x {n x ,t) 

For the simplest models with r x = const, we have 

d[n x (l + zf] 



dt 



£ x (t)H(t){n b ) , (5) 



1 f t-t u \ n x (t u ) 

where t u is the age of the Universe. 

Obviously, = /Cj£ x , e a = n a e x where factors and K a characterize the efficiency of 
resonance and ionized photons production. The decays of x-particles cannot significantly 
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distort the thermodynamics of the Universe at high redshifts z > 10 3 but this energy 
injection changes the kinetic of recombination at z ~ 10 3 that leads to the observable 
distortions of the power spectra of CMB anisotropy and polarization. To evaluate these 
distortions we firstly have to follow the transformation of high energy injected particles to 
UV photons influenced directly on the recombination process. 

The electromagnetic cascades initiated by the ultra-high energy cosmic rays 
(UHECR) are composed of photons, protons, electrons, positrons and neutrinos. 
At high redshifts, the cascades develop very rapidly via interaction with the CMB 
photons and pair creation ('Juhecr + Icmb e+ + e ~), proton-photon pair production 
ipuHECR+lCMB -> p'+7'+e++e~), inverse Compton scattering {e^HECR+ICMB -> e'+l), 
pair creation (e^} HECR + ^cmb — > e' + e~ + e + + 7'), and, for neutrino, electron- positron 
pair creation through the Z-resonance. As was shown by Berezinsky et al. (1990) and 
Protheroe et al. (1995), these processes result in the universal normalized spectrum of a 
cascade with a primary energy E 7 which can be written as follows: 



where E c ~ 4.6- 10 4 (l + ^)" 1 GeV, E a = 1.8-10 3 (l + z) -1 GeV. At the period of recombination 
z ~ 10 3 and for lower redshifts both energies, E a and E C) are larger than the limit of 
the electron-positron pair production E e + t£ - = 2m e = 1 MeV and the spectrum (7) 
describes both the energy distribution at E > E e + >e - and the injection of UV photons with 
E E e + y£ -. However, the spectrum of these UV photons is distorted due to the interaction 
of photons with the hydrogen - helium plasma. 

Let us briefly discuss the relationship between the initial spectra of the SHDM decay 




(7) 







E c < E 



v 
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and the parameters e a and Si of hydrogen excitation and ionization. Evidently, the energy 
density of UV photons is proportional to the rate of energy injection, Q(t), 

N(E) oc Q{t){E a /Ef' 2 = oo(t)Q(t u )(E a /E) 3 / 2 , (8) 



POO 

Q(t) = u(t)Q(t u ) = / dE EN(E) , 
Jo 



where Q(t u ) is the rate of energy injection at t — t u . For the net of ionizing UV photons 
with E > I ~ 13.6 eV we have 

Using normalization Q(t u )H~ l (t u ) ~ e EGRET to the EGRET energy density € E gret — 
4 ■ 1CT 7 eVcm~ 3 (Berton, Sigl&Silk 2001) we can estimate the functions Si and e a . Namely, 
for 

u(t) oct 4 ^oc (l + z)f (4 - p) (10) 

we get 

s.^e^ i(z) ( J , t J° B ! T 0) ) (1 + oc (1 + (11) 

The power index p — 1 corresponds to the model involving the release of x-particles 
from topological defects, such as ordinary cosmic strings, necklaces and magnetic monopoles 
(see, e.g., Sigl 2001). The model with p = 2 corresponds to the exponential decay of SHDM 
particles described by Eq. (3) with r x » ty. The model discussed in Peebles, Seager and 
Hu (2000) with = const, and e a (z) = const, corresponds to p — 4/3. 

For some models such as those of evaporation of primordial black holes or decays 
of SHDM relics with shorter life time, Tx = t(z x ) < t u , more complex normalization is 
required. Thus, for rx — 0.1t u that corresponds to redshift z x ~ 6 — 7, these decays occur 
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before or just at the epoch of galaxy formation and now their number is small. For such 
models and for z >> z x we get, instead of Eq. (11), 

Uod(z) = S(z)G(z x ) = £(z) exp ((1 + * x )§) . (12) 

For some models with different p and Q(z x ), functions e a (z) are plotted in Fig. 1. As 
one can see, for some of the models we should have e a > 10~ 3 — 10~ 2 at z < 1500 and 
distortions of the kinetics of hydrogen recombination can be quite strong. For other models 
these distortions are small and cannot be observed. 

Other approach to restrictions of admissible energy injection is discussed in Sec. 6 . 

2.2. Acceleration of recombination 

The basic idea of accelerated recombination is very transparent (Naselsky and Novikov 
2002). Let us consider the model in which a fraction of baryons is concentrated in high 
density low mass clouds. The hydrogen recombination inside the clouds proceeds faster 
than in the low density intercloud medium, and the mean ionization fraction is 

(x e \Z)) = x e ^ n Z m + a? e ,out(l Z m ) (13) 

where Z m is the mass fraction accumulated by clouds, x e ^ n and x e>ou t are the ionization 
fractions within and between clouds. This expression clearly shows that in such models the 
redshift variations of (x e (z)) differ from those in the standard models with homogeneous 
spatial distribution of baryons and at the same values of cosmological parameters. 

Cosmological models with such clumpy distribution of baryons appear when we consider 
the mixture of adiabatic and isocurvature (or isotemperature) perturbations (Dolgov & Silk 
1993). In suitable models isocurvature perturbations provide the concentration of some 
fraction of baryons within low massive clouds with M ~ 10 2 — 1O 5 M while adiabatic 
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perturbations dominate on larger scales and are responsible for galaxy formation and other 
effects. 

The ionization history of hydrogen in such models can also be described by Eqs. (1,2) 
with e a < and Si > (Naselsky and Novikov 2002). Of course, this choice does not mean 
that Lj-a photons disappear or are really transformed to Ly-c ones. But the suitable choice 
e a < allows to describe correctly the faster hydrogen recombination within the clouds 
while the choice Si > allows for the excess of Ly-a photons and the delay of hydrogen 
recombination in the low density intercloud medium. 

2.3. Three scenarios of recombination distortions 

For wide set of models, the functions e a and £j can be written as follows: 

^(^) = ^expK m 1cfO-, C = (1 + *)/(! + *d) , (14) 
j 

where parameters c- a , z&, rrij , and rij characterize the j th source of distortion. It is natural 
to expect that Si(z) > for all redshifts z < 10 3 . This additional ionization suppresses 
partly the CMB anisotropy because of the growth of the optical depth for the Compton 
scattering at z < 10 3 . 

Three scenarios can be discussed depending upon the shape of the function e a (z). 
Namely, the models with e a (z) > and e a (z) < describe the delay or acceleration 
of recombination, respectively. More complex models with the acceleration of the 
recombination at z > z ad and the delay of recombination at z < z ad are described by 
complex functions with e a (z) < at z > z a a and e a (z) > for z < z a( i- 

All these scenarios can be based on the realistic physical models. For example, for the 
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decay of the long lived SHDM particles we get 

e a ~ £j oc (1 + z)' 1 , (15) 

which corresponds to the first scenario. The same scenario is realized for the decay of 
topological defects and evaporation of primordial black holes discussed above. In all the 
cases, the Ly-c and Ly-a photons are formed via the electromagnetic cascades and e a ~ 
Third scenario is realized for the clumpy baryonic model discussed above. The list of the 
physical models can be essentially extended. 



3. Distortion of the CMB anisotropy power spectrum 

As was mentioned in Peebles, Seager and Hu (2001), the distortions of recombination 
shift the redshift of the last scattering surface, z r , and change its thickness, A z . In turn, 
these variations shift positions of the peaks in the CMB anisotropy and polarization power 
spectrum (Hu and Sugiyama 1995) and change their amplitudes. Here we give the rough 
analytical estimates of these influences. More accurate numerical results are presented 
below. 

Both parameters of the last scattering surface, z r and A z , are roughly expressed 
through the Thomson optical depth, tt, and its derivations, r' T = cLtt/cLz : 

TT -J rAz}dz > Tt[z) ~ h{z){i + z) -ITT' (16) 

Here ot is the Thomson cross-section and c is the speed of light. 

The redshift of the last scattering surface, z r , is defined by the position of the maximum 
of the so-called visibility function, g(r T ) = r^exp(— tt), 



dg(r T ) 
dz 



0, Tr(z r ) = (r' T {z r )} , (17) 
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and its width is 




(see, e.g., Naselsky & Polnarev 1987). Taking into account that recombination is a fast 
process and A z ~ 0.1z r , we can further consider only redshift variations of the ionization 
fraction, x e , and take the function $ in (16) at z = z r . This means that the parameters z* 
and A z (z*) can be found from the following equations: 

x' e (z*) ~ x 2 e (z*)$(z r )/z r , A z ~ x e (z*)/x' e (z*) m x~ 1 (z*)z r /^(z r ) . (19) 

Here z r and z* are the redshifts of the last scattering surface in the standard and the 
distorted models. For small distortions of ionization fraction we have 

X e ~ X e ^ s f -\- X e ^Z T {\ Zifj Z'f ) , X e ~ -^e,st dfl a / dz — %e,st ^a(^r)/ %r i 

z* ~ ^ r (l - e a (zr)i/), i/ = 0.5$~ 2 (z r )x-^ , (20) 

A z (z*) ~ A Zi8t [l + O^^)^ 1 ^)^] . 

Here x 6iSt and A ZjSi are the ionization fraction and the thickness of the last scattering 
surface for the standard model. 

The leading order dependence of the first peak location is (Peebles, Seager and Hu 
2001) 

Zi(e a ) ~ h,8t\/z*/z r ~ h iSt [l - 0.5ue a (z r )] . (21) 

For the second peak of the CMB power spectrum the shift of the position is practically the 
same as for the first one, but the ratio between their amplitudes is a linear function of ratio 
z*/z r . For the next few peaks growth of the thickness of the last scattering surface, A z , 
increases the damping and decreases the amplitudes of the peaks (White 2001). In contrast, 
for models with accelerated recombination and e a < 0, peaks are shifted to larger / and 
the damping of perturbations decreases together with A z as compared with the standard 
model. 
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4. Five models of recombination and distortions of the ionization history. 



To illustrate the possible manifestation of the external sources of resonance and ionized 
photons and non-uniform distribution of baryonic matter in the ionization history of the 
Universe we compare four models with different ionization history with the basic spatially 
flat ACDM model, Vl K = 0, the standard ionization history, Si = e a = 0, and 

tt A = 0.7, tt c + tt b = 0.3, h = 0.7, n s = l, Vt b h 2 = 0.022, Vt c h 2 = 0.125. (22) 

This baryonic matter density is close to the most probable value found from the comparison 
of the Big-Bang nucleosynthesis prediction with the observable abundance of the light 
chemical elements. The basic model includes also the optical depth for the Compton 
scattering due to reionization of the Universe by young galaxies and quasars at redshift z rei 



where x e ~ 1 is the fractional ionization at the redshift z re i and Y p is the mass fraction of 
helium. For parameters (22), Y p = 0.24, z re i ~ 14 and r re j ~ 0.1, this model provides the 
best-fit for the recent CBI data. We will refer to the basic model as to Model 1 . 

Recent observations of the farthest quasars (Fan et al. 2001) give the low limit for the 
epoch of reionization, z rei > 6, which does not contradict the adopted value z rei ~ 14. For 
z rei — 6, the Thompson optical depth is only r rei ~ 0.04. We will refer to the model with 
parameters (22), z re i ~ 6 and r re i ~ 0.04 as to the Model la and will discuss it in Section 6. 

To demonstrate the influence of distortions of the ionization history, we consider three 
models with various distortions. In Model 2 we change the baryonic fraction and accept 
Vl b h 2 = 0.03 preserving the total density, £l b + fi c = 0.3. The Model 3 is the cloudy 
baryonic model with {VL b )h 2 = 0.022 and reionization of hydrogen at redshift z rei ~ 14 



4.1. 



Models of recombination. 




(23) 
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as for the Model 1. We assume also that the density contrast between clouds, p in , and 
interclouds medium, p out , is rj = Pi n /p ou t = 11; clouds occupy the volume fraction /„ = 0.1 
and accumulate the mass fraction Z m = r]f v [rjf v + (1 — = 0.55. For simplicity, we 

consider all clouds with the same mass in range 1O 2 M < M d < 1O 5 M (Naselsky and 
Novikov, 2002). 

The Model 4 differs from the Model 1 by action of external sources of ionization with 

£ a = exp , e< = PC^ exp , 

Z = (l + z)/(l + z d )~10~ 3 (l + z). 

The choice a ~ 0.3, /3 ~ 0.13 coincides with the model of primordial black hole evaporation 
at the moment of recombination (Naselsky 1978; Naselsky & Polnarev 1987). In this Model 
we exclude the impact of late reionization which mimics roughly the model with z re i ~ 6. 
For some applications we consider also the Model 5 which differ from the Model 4 by the 
baryonic density, VL^h 2 = 0.03 preserving the total density, fi& + Jl c = 0.3. In this respect, 
the Model 5 is identical to the Model 2. 

Main parameters of the Models are listed in Table 1. 

For all the models the anisotropy and polarization power spectra are found with 
modified RECFAST and CMBFAST codes. To characterize the differences between the 
models and to compare them with the expected sensitivity of the PLANCK experiment we 
describe the deviations of the anisotropy and polarization power spectrum in terms of the 
functions 

DUi) = 2 [cm - q(0] / [cm + , 

Dim = 2 [cf (o - c m i i c ! (o + c m > ^ 

where i,j relate to the Models and C a and C p are the power spectrum of anisotropy and E 
component of polarization, respectively. The differences between the ionization fractions in 
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the models are characterized by the ratios 



Axu = ■2±—± . (25) 



4.2. Distortions of the ionization history 

The redshift variations of the relative difference between the ionization fractions in the 
basic Model 1 and three other models, Axij, j = 2,3,4 are plotted in Fig. 2 together with 
the redshift variations of the ionization fraction in the basic Model 1. 

For the Model 2, this difference does not exceed 20% at z ~ 10 3 and increases up to 
30% at lower redshifts which demonstrates a relatively weak influence of the variations of 
baryonic fraction. For the Model 3 with accelerated recombination, the relative difference 
Axi i3 < at redshifts z > 700 and the ionization fraction decreases with time faster than in 
the basic model. For z < 700 the ionization fraction is larger than that in the basic model, 
Ax 1)3 > 0, but the excess is less than 10% and the distortions of the CMB anisotropy and 
polarization power spectra remain small due to the small optical depth. 

The conspicuous delay of recombination and distortions of the ionization fraction 
appear for the Model 4. They achieve ~ 20 — 30% already at z ~ 10 3 and are maximal 
at z ~ 500, which provides the essential growth of the optical depth for the Thomson 
scattering and suppresses the amplitude of the CMB anisotropy. 

These results indicate that we can expect noticeable variations of the CMB anisotropy 
and polarization power spectra for the Model 3 and the Model 4 while for the Model 2 they 
can be quite moderate. 
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5. Anisotropy and polarization power spectra as a test of the ionization 

history of the Universe. 

To obtain the power spectra of the CMB anisotropy and polarization in our models we 
have to use the modification of the CMBFAST code (Seljak and Zaldarriaga, 1996) taking 
into account more complicated ionization history of the Universe discussed above. For the 
reference Model 1 and 2, we use the standard CMBFAST code for the optical depth r r and 
the ionization fraction x e — 1 achieved at the redshift of reionization z re i. 

In Fig. 3 we plot the CMB power spectrum for the models 1, 3 and 4 in comparison 
with the observational data from the BOOMERANG, MAXIMA-1, and CBI (Mason et 
al, 2002) experiments at the multipole range I < 2 000 where the possible influence of the 
Sunyaev-Zel'dovich clusters is not yet crucial. 

As one can see from Fig. 3, for all models the power spectra are very similar to each 
other and the Model 1, 3 and 4 are consistent with the available measurements of C(l). 
This means that the more complicated ionization history of hydrogen violates the standard 
fitting procedure. In particular, the comparison of results listed in Table 2 for the Model 
1,2, and 5 shows that the determination of baryonic density from the CMB measurements 
is ambiguous. This fact is important for the interpretation of the MAP and PLANCK 
measurements. 

To characterize the differences between the Models and to compare them with the 
expected sensitivity of the PLANCK experiment we plot in Fig. 4 the functions Dis(l) and 
-D14 (/) defined by (24) for the multipole range 2 < I < 2000. These functions can be directly 
compared with the expected error bars of C (I) for the PLANCK mission. 

As is clear from Fig. 4 for the multipole range 500 < / < 2000, 0.05 < D 13 (l) < 0.2 and 
0.05 < Di A (l) < 0.3, and there is a quite regular tilt at high multipole range. Obviously, the 
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w = (a p 9 FWHM ) 2 , Wi ~ exp 



fsky ^ 0.65 . 



tilt is related to the damping of the C(l) power spectrum in the basic model because of the 
late reionization with r rei = 0.1. In order to clarify the manifestations of the more complex 
ionization regimes in the Model 3 and 4 we have to compare the peak-to-peak variations of 
amplitudes, D 13 (l) and D u (l), with the expected errors of measurements for the PLANCK 
experiment. If the systematic effects will be successfully removed, the error should be close 
to 

^ffl - . 1 [1 + w^C-\l)Wr 2 } (26) 

1(1 + 1) 
212 

Here f s k y is the sky coverage during the first year of observations, a p is the sensitivity per 
resolution element Ofwhm x @fwhm and l s = \fS\n29p\ VHM . 

For example, for the HFI 353 GHz channel the #fwhm is 5 arcmin and the 
error at / < 2000 are defined by the first term in Eq (26) . As one can see from 
Fig. 4, at I ~ 1000 — 1500 the peak-to-peak amplitudes are ~ 10 — 20%, while 
AC(l)/C(l) ~ 5 — 7%. This means that such peculiarities of the anisotropy power spectrum 
can be detected by the PLANCK mission, if the systematic effects would be less than 
AC(l)/C(l)^[f sky (l + l/2)]-^. 

The polarization power spectra are even more sensitive to distortions of the ionization 
history of the Universe. Recently the DASI experiment (Kovac et al. 2002) has detected the 
CMB polarization without significant manifestation of polarized foregrounds. It indicates 
that the launched MAP and the PLANCK missions will be able to measure the polarization 
of the CMB with unprecendent accuracy. 

The polarization power spectrum AT p 2 (/) = 1(1 + l)C p (l)/2ir are plotted in Fig. 5 for the 
Model 1, 3 and 4 and their relative variations, D\-(l), with respect to the basic model are 
plotted in Fig. 6 . As one can see the differences between spectra are quite small especially 
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at I ^> 200. But the bump at I ~ 5 — 10 which is caused by the late reionization in the basic 
model can be essentially suppressed for the Model 4 wherein the required optical depth is 
gained at redshifts z ~ 100 — 500. 

For the PLANCK mission, even if all the systematic and foreground contaminations 
are removed from the polarization map (see Franco, Fosalba and Tauber 2002), one of the 
most important sources of uncertainty at I <C 200 is the cosmic variance. For f sky ~ 0.65 
and I ~ 2000 the cosmic variance limit is close to AC(l)/C(l) ~ 3% for la confidence level 
and all peculiarities of the polarization power spectrum should be observable. As it follows 
from Fig. 6, all polarization features which corresponds to the Model 1, 3 and 4 can be 
observed by the PLANCK mission, and, for the low multipole part I < 200, can be probably 
observed already by the MAP mission. This means that future polarization experiments 
can be crucial for investigation of the ionization history of the cosmic plasma as near the 
redshift of the last scattering surface as at the redshift of reionization. 

6. Low multipole peculiarities of the polarization power spectrum in models 

with delayed recombination 

The signature of late reionization is the bump at small multipoles, I < 10, in the 
polarization power spectrum. Its amplitude and shape depend upon the ionization history 
characterized by z rei and r rei . This bump is clearly seen in Fig. 5 for the basic model with 
z rei = 14 and r rei = 0.1 but it is strongly suppressed for the later reionization such as, 
for example, for z re i ~ 6 and r re j = 0.04 widely discussed in connection with the recent 
observations of the farthest quasars (Fan et al. 2001). For the non-standard ionization 
history discussed in this paper, the shape of the polarization power spectrum is more 
complicated but its observations allow to restore some features of the ionization hystory. 
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As is seen from Fig. 6, for the Model 3 with accelerated recombination this bump is 
just the same as in the basic Model I. However, in the Model 4 with delayed recombination 
and without any sudden reionization at low redshifts significant variations of D^j(l) are 
seen in the Fig. 6 at I < 500. They can be observed by the PLANCK experiment even for 
relatively small Ei(z), which would clarify the ionization history for such models. 

To illustrate these abilities, we plot in Fig. 7 the relative distortions of the polarization 
power spectrum at I < 1000 for several models with delayed recombination and different 
£i(z). First of all, to illustrate the influence of the optical depth, r re j, we compare the Model 
1 with the Model la with the same parameters and z rei ~ 6 and r rei = 0.04. Large values 
Df j(l) for these models at / < 20 allows to easily discriminate between these models. 

In the same Fig. 7 we compare also the Model la with several models of delayed 
recombination with the same cosmological parameters (22) and with 



where z r ~ 10 3 is the redshift of the last scattering surface. We see that for the models 
with p — 1, Szei(z r ) = 10~ 2 and with p = 2, hei{z r ) = 10~ 3 the relative distortions, 



p = 1, &,Ei(z r ) = 10~ 3 distortions become small and are under observational limits. These 
comparisons demonstrate the abilities of feature polarization measurements with respect to 
restrictions of the ionization history of the Universe at z < 10 3 . 

In wide class of such models the injection of ionized photons at z < z r increases the 
ionization fraction and the Thompson optical depth with respect to the standard models. 
In such case, the already available limitations of the depth t(z) < r t h r ~ 0.1 restrict the 
admissible £i(z) even at high redshifts where estimates (11) are not effective. To estimate 
the optical depth for such models we have to trace numerically the ionization fraction at 
all redshifts. However, even rough analytical estimates are informative and allow to reject 




(27) 



Dij(l), are quite similar and can be observed at / < 100. However, for the model with 



-20- 



some models. 

For such rough estimates we use the equation of ionization balance 

a c (n b )x 2 e ~ Si{z)H{z) (28) 

where a c oc T~ 0M is the coefficient of hydrogen recombination and T m is the temperature 
of the matter. As is well known, T m oc 1 + z for z > z eq ~ 10 2 and T m oc (1 + z) 2 for lower 
redshifts and, so, a c oc T~ - 6 oc (1 + z)~ 2uJ , oo ps 0.3 for z > z eq and 0.6 for z < z eq . 
Thus, for Si{z) given by (27) we get for the ionization fraction and for the optical depth 

x e (z) oc ef(z r )(l + zY, n = -0.75 +u + (1 - 3p/4) , 

,,..36 ri±±y (i^) mv 2 f^iy 1 ' 4 (2 ») 

1 v rJ \i + Zr) V o.76 yvo.o2y Vo.125 y v ; 

where Y p is the mass fraction of helium, v — 1.5 + /i — 0.75 + + (1 — 3p/4), and z T is the 
maximal redshift where the ionization balance is described by (28). This relation is directly 
applied for z T < z eq , oj = 0.6 and can be easily extended for the case z T > z eq . 

The observations of farthest quasars show that the reionization occurs at z rei > 6 when 
r ~ 0.04 and, so, for cosmological parameters (22), Y p = 0.24 and r thr ~ 0.1 we get 

e t (z r ) < 10''u 2 (j^j ^ ■ (30) 

As compared with the normalization (11), this expression more effectively restricts the 
possible energy injection at high redshifts. 

7. Discussion and conclusion 

The already achieved precision of measurements of the CMB anisotropy and especially 
expected one for the MAP and PLANCK missions allows to discriminate wide set of 
cosmological models and to reveal noticeable distortions of the standard model of hydrogen 
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recombination process at redshifts z ~ 10 3 and of the ionization history. Many physical 
models can be considered as a basis for such discussions and can be tested and restricted 
via such measurements. 

In this paper we compare available observations of the CMB anisotropy with those 
expected in several models with different kinetics of recombination. The differences can be 
caused by possible external sources of resonance and ionized radiation, possible small scale 
clustering of baryonic component and other factors. We show that the models with more 
complicate ionization history provide better description of the available observational data. 
This means also that the accepted estimates of cosmological parameters can be changed 
when such models are incorporated in the fitting procedure. It is especially crucial for the 
estimates of baryonic density, fl^h 2 , which is closely related to the Big-Bang nucleosynthesis 
predictions. 

In this paper we do not consider the impact of Sunyaev-Zel'dovich (SZ) effect and the 
weak lensing by the large scale structure elements. However, the SZ effect is frequency 
dependent that allows to reduce its contribution for multifrequency measurements. As was 
shown in Schmalzing, Takada and Futamase (2000), the influence of lensing effect can really 
distort the measured characteristics of the CMB only at larger I > 10 3 while the discussed 
deviations of ionizing history distort these characteristics at all /. This means that the 
measurements of the CMB spectra in wide range of / allows to discriminate between these 
distortions. However, this problem requires more detailed discussions for various models of 
both lensing and recombination. 

Near the first Doppler peak, the cosmic variance effect restricts an accuracy of estimates 
of C a (l) as AC a (l)/C a (l) ~ 10% for all experiments with the sky coverage f sky ~ 65%. 
For the recent balloon-borne and ground-based experiments the achieved accuracy in this 
range of / is even worse than this limit. With such precision, we cannot discriminate the 
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basic model and models with delayed recombination and another Vt b h 2 . To illustrate this 
statement we compare quantitatively the observational data with the power spectra for 
Model 1, 2&3 parameters which are given in Sec. 4.1 . Instead of the Model 4 we here 
present the results for Model 5 which differs from Model 4 by the values of parameters 
VL b h 2 = 0.03 and VL m = 0.256 which are the same as in the Model 2. 

The quality of models in comparison with the recent BOOMERANG, MAXIMA-1, 
CBI and VSA observational data can be characterized quantitatively by the x 2 -parameter 
listed in Table 2 for all models. Here the reference Model 1 provides the best- fit for the 
CBI data for the standard ionization history while xi > x\ f° r the Model 2 with larger 
VL b h 2 . However, the Model 3 demonstrates an excellent agreement with the CBI data, 
x\ = 0.5xi, and xi < Xi f° r other observations. For the Model 5 with Vt b h 2 = 0.03 we get 
also xi < Xi for the CBI data and only for the BOOMERANG data xi > xl 

Obviously, such extension of the cosmological models increases the number of 
parameters using to fit observed power spectra of the anisotropy and, in perspective, of the 
polarization. We show that the influence of these "missing" parameters, namely Si{z) and 
e a (z) , can improve the fits used and, in particular, to obtain better agreement between the 
observed and expected positions and amplitudes of peaks for higher multipoles. At the same 
time, including of the "missing" parameters changes the measured values of cosmological 
parameters usually discussed. 

We show also that the expected sensitivity of the MAP and the PLANK missions in 
respect of the measurements of the polarization will allow to discriminate between main 
families of such models and, in particular, between models with small and large optical 
depth at redshifts 10 3 > z > 20-50. We would like to note that the realistic values of the 
cosmological parameters can not be obtained from the CMB data without the PLANCK 
observations of the polarization. 
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Fig. 1. — Redshift variations of e a for models (11, 12) for p — 0, 0(;z x ) = 1 (solid line), 
p = 2, e(z x ) = 1 (dot-dashed line), p = 2, Q(z) = 1 (dot -dashed line), p = 2, 9(z) = 10 6 
(dot-dot-dot dashed line) and p = 2, 6(z) = 3 • 10 5 (dashed line). 
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Fig. 2. — Redshift variations of relative differences of the ionization fractions between basic 
Model 1 and other Models, Axij, for j = 2,3&4 (thin solid, dashed and dot-dashed lines, 
respectively). Thick solid line shows the variations of the ionization fraction for the basic 
Model I. 



-30- 




Fig. 3. — CMB anisotropy power spectrum for the Model 1 (solid line), the Model 3 (long 
dashed line) and the Model 4 (dot-dashed line) in comparison with the BOOMERANG 
(points), MAXIMA -1 (large X) and CMI Ml (small x) data. 




Fig. 4. — Relative differences of CMB polarization power spectra, D® 3 (1) (solid line) and 
D" 4 (l) (dashed line). Thick solid line show expected errors for the PLANCK mission (with- 
out systematics errors). 
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Fig. 5. — CMB polarization power spectra for the Model 1 (solid line), the Model 3 (dashed 
line), and the Model 4 (dot-dashed line). 
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Fig. 6. — Relative differences of the CMB polarization power spectrum, D\ -, between the 
basic Model 1 and the Model 3 (dashed line) and the Model 4 (solid line). Thick solid line 
show expected errors for the PLANCK mission (without systematics errors). 



-34- 




Fig. 7. — Relative differences D\ k of the CMB polarization power spectra for low multipoles 
for four different models of ionization. The difference between the basic Model 1 with 
Trei = 0.1 and Model la with the same parameters and r re j = 0.04 is plotted by solid 
line; the differences between the Model la and models (27) with p — 1, e(z r ) = 10~ 2 and 
e(z r ) = 10~ 3 are plotted by dashed and dot dashed lines; the difference between the Model 
la and model (27) with p — 2, e(z r ) = 10~ 3 is plotted by dot line. Thick solid lines show 
expected errors for the PLANCK mission. 
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Table 1: Main parameters of the Models 









n b h 2 




Trei 


1 


basic 


0.256 


0.022 


14 


0.1 


la 


basic 


0.256 


0.022 


6 


0.04 


2 


basic 


0.240 


0.03 


14 


0.1 


3 


accelerate 


0.256 


0.022 


14 


0.1 


4 


delayed 


0.256 


0.022 





0.0 


5 


delayed 


0.240 


0.03 





0.0 
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Table 2: x 2 for five observed anisotropy power spectra and for the Models 1, 2, 3, &5 





BOOM 


MAX 


CBIM1 


CBIM2 


VSA 


1 


11.2 


14.1 


3.20 


7.15 


7.1 


2 


34.6 


17.7 


2.75 


7.51 


13.0 


3 


8.7 


11.3 


1.60 


5.62 


6.6 


5 


27.7 


15.7 


1.87 


5.15 


10.0 



